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Abstract 

The Lagrangian and Hamiltonian formulations of electromagnet ism are reviewed and the 
Maxwell equations are obtained from the Hamiltonian for a system of many electric charges. It is 
shown that three of the equations which were obtained from the Hamiltonian, namely the Lorentz 
force law and two Maxwell equations, can be obtained as well from a set of postulated Poisson 
brackets. It is shown how the results derived from these brackets can be used to reconstruct the 
original Lagrangian for the theory aided by some reasoning based on physical concepts. 

PACS: 45.20Jj, 45.10Db, 02.30. Zz, 03.50.De 



1. Introduction. 



The Maxwell equations provide a type of mathematical summary of several fundamental laws of 
electromagnetism which originally had their origin in experimental observations. The subsequent 
development of a non-Abelian extension of the Maxwell theory to the non-Abelian Yang-Mills 
form and the diverse applications of Yang-Mills to the subatomic realm of high energy physics has 
probably given further impetus to the study of Maxwell's theory. 

A very novel derivation of a pair of the four Maxwell equations was originally introduced by 
Feynman, but the exact details of his argument remained unpublished until some of the essential 
arguments were presented by Dyson [1]. It has also been shown that this procedure can be 
generalized as well to the case of the dynamics of particles which possess other internal degrees 
of freedom [2]. The idea at the heart of these derivations is to postulate a fundamental set of 
Poisson brackets between the fundamental variables of the system. The basic defining relations 
of the bracket are then applied to the original collection of brackets as well as other operations, 
such as differentiation, to generate further new relations and connections between the variables of 
the problem [3] . This process is capable of generating some of the Maxwell equations, as we will 
show. 

In this paper, we begin by developing the Hamiltonian formulation of the Maxwell theory, and 
derive the Maxwell equations in this context. It is shown that the Hamiltonian formalism of the 
classical system, like the Lagrangian formalism on which it is based, is also invariant under gauge 
transformations [4-5] . Different Hamiltonians can be written so that they all have the same form, 
often referred to as minimal electromagnetic coupling. Next, the development of the Maxwell 
equations from a minimal set of defining brackets involving the variables and the basic algebraic 
properties of the Poisson bracket is reviewed [3]. To this end, a fundamental set of brackets is 
postulated at the outset in a natural way, and subsequently the algebraic properties of the bracket 
as well defined analytic operations are applied to obtain the basic Lorentz force law as well as a pair 
of the Maxwell equations. Based on these initial results, some of the ideas of the inverse problem 
of the calculus of variations are applied [6] . To emphasize, the procedure relies on defining a basic 
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set of brackets and using fundamental properties of the bracket to generate new relations, such as 
the Leibnitz rule and Jacobi identity, regardless of the underlying definition of the bracket It is 
shown from these results and with the help of some additional physical motivation at the end that 
the full Lagrangian for the theory can be reconstructed. In the sense that a Lagrangian theory 
can be formulated out of a set of elementary results, the complete set of Maxwell's equations can 
be obtained. The results that are obtained from the Lagrangian by means of the Euler-Lagrange 
equations, can then be used to define a Hamiltonian for the theory to complete the construction 
[7-8]. 

It may be asked why this approach is adopted. There are several approaches already that begin 
this kind of development with commutators [1]. Here we show it is possible to proceed entirely 
in the classical domain. Of course, the Maxwell equations exited before and independently of 
quantum mechanics, and nonetheless they are fundamental in generating quantum theories of 
electromagnetism. The same type of analysis can be carried out on classical theories of gravity 
and it may prove possible to adopt some of the ideas here to proceed to quantum theories of 
gravity. 

2. Lagrangian and Hamiltonian Formalism. 

The system which is of interest here consists of a collection of nonrelativistic particles which 
interact with an external electromagnetic field. The Lagrangian for the system is sufficient to be 
used with the principle of least action to generate the equations of motion. Moreover, a Lagrangian 
is required to construct a Hamiltonian for the system. 

A system of nonrelativistic particles, each having a charge e a , mass m a and a displacement 
vector r a (t) at time t for a — 1, 2, • • • , N in an electric field E(x, t) and magnetic field B(x, t) can 
be described by the Lagrangian L, which in Lorentz-Heaviside units takes the form 

L= l -J ^(E 2 (x,t)-B 2 (x,t)) + if;m a r2+i J d 3 x (J • A - cpA ) - U. (2.1) 

For the sake of generality, an arbitrary static external potential energy U can be included in the 
Lagrangian as well, but it is not essential for what follows. The charge density p(x, t) of the system 
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(2.1) is defined by 



p(x,*) = ^e Q <J(x-r a (t)), (2.2) 



a=l 



and the current density J(x, £) associated with the motion of the particle is given by 



JV 



J(x,t) = ^e a r a 5(x-r Q (t)). (2.3) 



a=l 



Of course, the charge density and current density satisfy the equation of continuity. 

The total electromagnetic field is characterized by the vector potential A and a scalar potential 
A Q in an arbitrary gauge. The total electric and magnetic fields are related to the potentials in 
the following way 

1 (9 A 

E(r, t) = -VA — , B = V x A, (2.4) 
c at 

where V is a vector operator and V a is the vector operator corresponding to particle a. The 
equations of motion can be obtained from Hamilton's principle of least action by varying the 
action S 



\ dtL. 
Jti 



The resulting Euler-Lagrange equations take the general form 

d dL dL ^ dL 

dtdQ, dQ j ^ l d(d t Q 3 Y 1 • ; 

where the Qi pertain to any of the physical variables on which L depends. For example, taking 
Q to be A , variation of the action S with respect to A gives Gauss's law. This can be obtained 
as well from the Euler-Lagrange system (2.5) by identifying Qi with A and determining the 
derivatives in (2.5). 

The Hamiltonian formalism for the total system is manifestly gauge invariant, and can be 
determined by first calculating the canonical momenta. The canonical momentum conjugate to 
the coordinate r a is given using (2.5) as 

Pa = T77- = mr a + —A(r a , t), (2.6) 

OT a C 

The canonical momentum conjugate to the field A is 



Since the Lagrangian is independent of the quantity A , the canonical momentum conjugate to 
A is n = 0. The Hamiltonian for the total system is defined by 

N 



h= [ d 3 x (n • A + n i ) + p« • *<* - L - 

J a=l 



(2.8) 



Substituting the canonical momenta II, Ilo and r a in terms of p a from (2.6), as well as the 
Lagrangian L given in (2.1), we obtain that 

N 

H = J2—P a -(p a --A)-l [ d 3 x (E 2 -B 2 ) -- [ d 3 x(E-A) 
m a c z J c J 



~\ E — (P« " - A ) 2 + U -- f d 3 x ( J • A - cpA ). 
2^m a c c J 

a=l 

Replacing A in the third term of (2.8) using (2.4), we can write 



(2.9) 



H = \( d 3 x(E 2 + B 2 ) + Ti(p a -^Af + [/+ ( d 3 xE-VA 
2 J 2 ^— \ m a c J 

N 

+ y ^ A • (p Q - ^A) - - / d 3 x (J • A - cpA ). (2.10) 
^— f m a c c c J 

Grouping the last three terms in (2.10) together, this equation can be rewritten using the expres- 
sions for p and J given in (2.2) and (2.3) 

i r N i r 

H=- d 3 x(E 2 (x,t) + B 2 (x,t)) + J2^(P«- e -^Mr«,t)) 2 + U+ / d 3 x(p-V-E)A . (2.11) 

A total derivative term has been dropped to write the Hamiltonian H in the form (2.11). 

From the Hamiltonian (2.11), Hamilton's equations can be developed. The first of Hamilton's 
equations for particle a is found by differentiating H with respect to p a , 

1 . c a , 

r a = - — = — (p a A. 2.12 

dp a m a c 

Using (2.12) and the vector identity V a (r a • r a ) = 2(r Q • V a )r Q + 2r a x (V Q x r Q ), we obtain the 
second of Hamilton's equations, 

OH e e dU 

p a = -— = — -i a ■ V a A(r Q , t) - —i a x (V a x A(r a , t)) + ^— + e a V a A (r a , t). (2.13) 

OT a C C OT a 



Solving for p a in (2.12) and differentiating this with respect to time, we find 

p a = m a r a + — (A a + r Q -^^). (2.14) 

C OT a 

Substituting (2.14) into (2.13), and simplifying with equation (2.4) yields 

e dU 
mr a = e Q E(r Q , t) + — r Q x B(r a , t) - j—. (2.15) 

C OY a 

This is Newton's second law in terms of the Lorentz force. 

It remains to complete this process for the field variables and the procedure is the same as for 
a system with a finite number of degrees of freedom. The functional derivative of H is calculated 
with respect to II defined by (2.7), and this gives the first of Hamilton's equations 

A = |g = -c(E + VA ), (2.16) 

where the time derivative of A is a partial derivative. 
The second of Hamilton's equations is given by 

n = "|x = - VxB + ^ (2 - 17 > 

which is the Ampere-Maxwell law, 

VxB = -J + -^. (2.18) 
c c at 

The equation for n is given by 

rio = -^ = -p + V-E. (2.19) 

The scalar potential in E given in (2.4) is not varied since E is proportional to II, hence indepen- 
dent. Now, since n = 0, equation (2.19) implies Gauss's law, 

V • E = p. (2.20) 

Finally, the equation for A is found by differentiating H with respect to n . Since n = 0, this 
equation is meaningless and does not provide an equation. This procedure has in fact generated 



all four of Maxwell's equations, although two have not yet been explicitly written down. To obtain 
the remaining two equations, first take the curl of (2.16) and use the identity V x VA = and 
(2-4), 

1 f) 1 r)B 

V x E = —V x VA — V x A = — . (2.21) 

cat c at 

This is Faraday's law. The condition that magnetic monopoles do not exist follows by taking the 
divergence of B given in (2.4), 

V • B = 0. (2.22) 
To summarize, the Maxwell's equations are given by (2.18), (2.20), (2.21) and (2.22). 

3. Maxwell Equations and Poisson Brackets. 

The algebra of classical observables on the manifold M will be denoted by T . A Poisson 
structure on a manifold M is a skew-symmetric bilinear map which is denoted {, } : T x T — > T 
such that (i) (JF, {, }) satisfies the Jacobi identity {G, {H, K}} + {H, {K, G}} + {K, {G, H}} = 0, 
(ii) the map X G = {, G} is a derivation of the associative algebra F(M) on M [9]. It satisfies the 
Leibnitz rule {G, H K} = H{G, K} + {G, H}K. A manifold M, which is endowed with a Poisson 
bracket on F(M), is called a Poisson manifold. These basic algebraic properties are used on their 
own without reference to a specific form for the bracket to develop the stated results. 

Let the local coordinate variables on the manifold be written in the form (w a ) = (q l , v l ), where 
i = 1,2,3. Indices will be raised and lowered in a trivial way using 5^, and repeated indices are 
summed over. Here, the q l may be interpreted as position coordinates and v % represent velocity 
components. The Poisson brackets are postulated in the following way [2-3], 

{qi,qj} = 0, m{q i ,v j } = 6 ij . (3.1) 

Any function H G T will define a dynamical system on M by the equation 

f = {G,H}. (3.2) 

Let us postulate an H e T such that equations of motion can be obtained from (3.2) as follows, 

q* = {q\ H} = v\ mv* = m{v\ H} = F\ (3.3) 



Differentiating the second bracket in (3.1) with respect to time generates the equation 

{QuVj} + {qi,Vj} = 0. (3.4) 
Multiplying both sides by m and substituting the equations of motion, there results the expression 

™{ft, ft} + {ft, F j} = 0. (3.5) 
Since the bracket is bilinear, this equation can be put into the form 

{{ft, Fj}, Ik] + m{{qi, qj}, q k } = 0. (3.6) 
Substituting q i: <jj and q k into the Jacobi identity, we have 

{{ft, Qj}, Qk} + {{<lj, Qk}, ft} + {{Qk, ft}, ft} = 0. (3.7) 
The bracket {ft,fto} is proportional to Sjk by (3.1), so (3.7) reduces to the constraint 

{{ft, ft}, ft} = 0. (3.8) 
Substituting (3.8) into (3.6), we obtain that 

{ft,{ft,^}} = 0. (3.9) 

The tensor {% Fj} is therefore antisymmetric due to the bracket property. This can be expressed 
in its dual form by the relation 

{ qi ,Fj} = -— e ijk B k (r,t). (3.10) 
mc 

Substituting (3.10) into (3.9), a bracket which contains qi and B k can be obtained 

{qi,B k } = 0. (3.11) 

The postulated relations (3.1) imply that the vector B depends only on the position and time of 
the particle. Equations (3.11) and (3.1) imply that Fj is at most linear in the velocities, and so 
we may write 

Fi(r,t) = eEi(r,t) + ^e ijk v j B k (r,t). (3.12) 



This is the Lorentz force law and serves to define the electric field. Using the property of bilinearity 
and the derivation property, we have 

6 6 

{ft, eEj + -e jak v a B k } = e{q i7 Ej} + -e jak {q h v a B k } 

6 6 

= e{q h Ej} + -e jak {qi,v a }B k + -e jak v a {qi, B k } (3.13) 

6 6 

= e{qi, Ej} H eji k B k = e{qi, Ej} e ijk B k . 

mc mc 

Since the result of (3.13) is the left-hand side of (3.10), it follows that 

{q t ,E 3 } = 0. (3.14) 

This bracket implies that the vector E, as in the case of B, depends only on the time and position 
coordinates. Equations (3.5) and (3.10) can be combined and will lead to a new equation for B k 
in terms of the bracket 

IT = ^c*{q it q j }. (3.15) 

Applying the Jacobi identity to the variables qi, <jj and q k and then contracting with e l: * k , there 
results 

^ k {qiAqj,qk}} = o. 

Replacing the bracket {qj,q k } in this using (3.15) generates a new bracket involving B, 

{qi,e ljk { qj ,q k }} = ^{q h B l } = 0. 

77TC 

This implies that {qi,Bi} =0, and therefore 

V-B = 0. (3.16) 

To obtain a second equation, let us begin with the equation for B s , given in (3.15). Differentiating 
both sides with respect to t, we have 

dB s dB s m 2 c sij . m2 c siU . ... m 2 c . 

ST + flfc* = 17 £ + 17 e lq " q ' } = T te ' 9} ' (X17) 
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Substituting (3.12) into the right hand side of (3.17), dividing out the common factor of e, and 
using property (ii) of the Poisson bracket, there results 

I ■ ■ TTl 

me SiJ {Ei + -e ial q a B h q t } = me stJ {Ei, qj } + — {B s , qj}qj 

771 771 771 

+ -{q,A 3 }B s - -W-'lAlh - -<Y{I[r<h\- (3-18) 

The second to last term on the right hand side of (3.18) is zero by symmetry, and upon substitution 
of the equation {q h B t } = into (3.18), this expression reduces to the following 

dB, dB, ., dE, ldB, 



s 



cdt c dqj €sjl dqj c dqj 

Simplifying this, the following Maxwell equation is obtained in the usual form 

1 <9B 

-— + VxE = 0. (3.19) 
c at 

4. Reconstruction of the Lagrangian. 

The fundamental Poisson brackets have been assumed to have the structure given in (3.1), in 
particular, if we combine (3.1) and (3.3), the basic relation is of the form 

m{q\q j } = 5 ij . (4.1) 

Consider the classical equations of motion with all masses set to unity of the form 

ft = /*(?, q,t). 

A nonsingular matrix Wy and a function L(q, q,t) are sought such that 

The conditions for the existence of Wij and L are called the Helmholtz conditions [10]. If a 
Lagrangian L exists, then is given by 

w « = wk- (43) 
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From (4.1), we take Wij to be proportional to 5ij. For the Hessian (4.3) to be invertible, the 
Lagrangian must obey 

mSij. (4.4) 



dtfdqi 

Integrating (4.4), it can be seen that if L exists, it must have the form 

1 e 1 e 

L = -mq % q 3 6 ii + -<f A { - eA + C = -mr 2 + -r • A - eA + C(E, B, A, A ). (4.5) 

2 c 2 c 

where the quantity C can be regarded as a constant of integration which will be a functional of the 
fields which were introduced in producing the solution (3.12) and equations (3.16), (3.19). Such a 
term would represent an energy contribution associated with these fields, and its exact structure 
can be determined next on physical grounds. The existence of L, however, does follow from the 
Helmholtz equations, as will be seen. If some physical ideas are now introduced and applied, the 
Lagrangian (4.5) can be generalized to the structure given in (2.1) by using (4.5) and writing 

N 



(4.6) 



L = J2 \ m ^l + ^fd 3 x(J-A- c P A ) + C(E, B, A, A ), 

where p and J have been defined in (2.2) and (2.3). 

We say that the force F(t, q> , q j , q>) is potential if Lagrange's equations of motion 

dT _ d_dT _ 

dq i dt dq 1 

are variational, that is, if there exists a Lagrange function L such that 

dT _ d_dT_ p _ dL_ _ d_dL_ 

dq i dtdq 1 1 dq l dtdq 1 

The necessary and sufficient conditions for a force F l (t, g J , g J ) to be potential can be written 

m + dF l = Q ^ 

dqi dq 1 

(4.7) 

dqi dq l dt dq 1 
The first of these can be differentiated with respect to qi to give 

d 2 F, d 2 F k d 2 Fj d 2 F t 



dqidq k dcfdqi dq i dq k dqidq k ' 
11 



therefore, 

° 1 = 0. 
dqidq k 

This equation can now be integrated to give the result Ft = a^q^+bi, where and b{ are functions 



of (t, g fc ). Substituting Fj into the pair (4.7), we obtain a set of three conditions on the and 6; 
as follows 

da is da sj d aji 

By setting a^- = —(e/c)eijkB k and eF l = 5 lj bj, we obtain the Lorentz force law (3.12) after 
reversing the direction of the velocity vector, such that E and B satisfy, 

ldB 

V-B = 0, VxE = — . 

c at 

These results are consistent with the development from the Poisson bracket point of view. It is 
interesting to note that if we make the opposite selection fej = S^B^ and = eijkE k , albeit for 
which there is no known force law, the remaining two sourceless Maxwell equations appear. Of 
course, the Maxwell equations are known to be symmetric under the transformations E — * B and 
B -> -E. 

To determine the quantity C in (4.6), we have to be allowed to reason from a physical point 
of view. This functional will be a scalar formed from the fields E and B, which will serve as a 
Lagrangian for these fields. To this end, we construct a functional of the form 

L em = y"d 3 x(aE 2 (x,t)+/3B 2 (x,t)+ 7 E(x,t)-B(x,t)), (4.9) 

and identify L em with the integration constant C. This reflects the fact that the fields involved 
will have an energy density associated with them in the absence of particles. Adding terms such 
as A 2 and A$ to this integrand would destroy the form invariance of the total Lagrangian under 
the addition of a total time derivative. 
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Since E is a vector and B is a pseudovector, this will be invariant with respect to parity only if 
7 = 0. The constants a and f3 can be taken so that other equations and quantities obtained from 
L will have their standard forms, in this case, a — (5 — 1/2. Now, from the Lagrangian and the 
principle of least action, we can immediately determine the remaining two Maxwell equations. If 
the action is varied with respect to A , Gauss's law results and if we vary with respect to A, then 
the Ampere-Maxwell law results as before. 

To summarize, it has been shown that a set of basic Poisson brackets leads to some of the 
basic structures in electromagnetism, in particular, the Lorentz force law. This can then be used 
in turn to develop a Lagrangian which must exist on account of the Helmholtz conditions. It is 
also quite interesting that without any direct appeal to Lorentz invariance, a system of equations 
results which can be shown to be relativistically invariant under Lorentz transformations [11-12]. 

It can be seen that the Maxwell equations are partitioned into two groups, and this is clearly 
indicated in the details of the reconstruction. Both Faraday's law (2.21) and (2.22) are kinematical 
laws. They follow from the relationship between the actual fields and the charged particles. It is 
not surprising then that these two equations are generated by the set of Poisson brackets [13]. On 
the other hand, the Ampere-Maxwell law (2.18) and Gauss's law (2.20) are dynamical equations. 
It would be of interest to know whether other types of equations of physical interest, for example 
to the realm of gravity [14], can be developed along similar lines as described here. 
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